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q . Abstract 

We reduce the perturbative uncertainty in the determination of \ V u b\ from inclusive 
semileptonic B decays by calculating the rate of B — > X u I v events with dilepton 
(D | invariant mass \fq* > Mb — Md at subleading order in the hybrid expansion, and 

• to next-to-leading order in renormalization-group improved perturbation theory. 

. We also resum logarithmic corrections to the leading power-suppressed contribu- 

tions. Studying the effect of different 6-quark mass definitions we find that the 
5_i ■ branching ratio after the cut is Br(5 — > X u lv) = (20.9 ± 4.0) 1 | 2 , where the 

dominant error is due to the uncertainty in the 6-quark mass. This implies that 
\V u b\ can be determined with a precision of about 10%. 



1 Introduction 



The total inclusive, semileptonic B — > X u I v decay rate admits a clean theoretical analy- 
sis based on the heavy-quark expansion. The optical theorem is used to relate the decay 
rate to the forward matrix element of the product of two hadronic currents, which is 
evaluated using the operator product expansion |l|, @, |J. The total rate is presently 
known including corrections of order a 2 M and A 2 /m\ (with A a typical hadronic scale). 
An accurate knowledge of the total rate does not immediately translate into a precise de- 
termination of \ V u b\, since experimentally one has to impose cuts on the kinematical vari- 
ables in order to discriminate the b — > u transitions against the much larger B —>■ X c l v 
background. The first inclusive determination of \V u b\ was based on a measurement of 
the charged-lepton energy spectrum in the endpoint region, which is kinematically for- 
bidden for decays with a charm meson in the final state. This cut eliminates about 90% 
of the signal. It is difficult to reliably calculate the remaining event fraction, since the 
operator product expansion breaks down in this part of the phase space and must be 
replaced by a twist expansion 0, ||. 

Recently, it was shown that the situation is improved when a cut q 2 > q$ with 
Qo — {Mb ~ M D ) 2 is applied to the dilepton invariant mass |7J]. This keeps about 
20% of the signal events, and the fraction of the remaining events can be calculated 
in a model-independent way using the operator product expansion. After the cut, the 
typical momentum flow fi c through the hadronic tensor is of order the charm-quark mass 
or less, so that the problem at hand involves the three scales mj > fi c > A. In the 
standard heavy-quark expansion, the charm and the 6-quark masses are treated on the 
same footing, but it is known that the convergence of this expansion at the charm scale 
can be poor. In this letter we follow the strategy proposed in || and perform a two-step 
expansion (the "hybrid expansion") in the ratios [i c / m b an d A//i c . The first step is to 
integrate out the degrees of freedom above the scale /i c and to sum up logarithms of the 
form a s ln(// c / nib) using the renormalization group (RG). This is achieved by expanding 
the weak currents and the QCD Lagrangian to order I /nib in the heavy-quark effective 
theory (HQET) ||, and calculating the Wilson coefficients of the pertinent operators at 
next-to-leading order (NLO) in RG-improved perturbation theory. Whereas the Wilson 
coefficients of the leading-order currents [K], |ll|] and of the 0(1 /nib) corrections to the 
HQET Lagrangian [12, [13| were known for a long time, the NLO calculation of the 



coefficients entering the expansion of the currents at order 1/rrib was completed only 



recently [14, 15, 16, 17 



In the present work, we perform the second step in the hybrid expansion and construct 
the operator product expansion in the resulting low-energy effective theory. We calculate 
the decay rate T(B — > X u I v)\ q 2 >q 2 at subleading order in the ratio ix c /mb and to NLO 
in RG-improved perturbation theory. We also resum leading logarithmic corrections 
to the power corrections of order (A//x c ) 2 . In contrast with previous analyses J7|, we 
calculate the decay rate directly, rather than calculating the fraction of events that passes 
the g 2 -cut and then multiplying this fraction with an independent result for the total 
semileptonic rate. Since the characteristic scale of the total rate is significantly higher 
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than (and parametrically different from) the typical hadronic scale ji c after the g 2 -cut, 
the theoretical error cannot be reduced by separating the two calculations. Based on a 
detailed analysis of the residual renormalization-scale dependence and the sensitivity to 
the definition of the 6-quark mass, we conclude that the B — > X u I v branching ratio after 
the dilepton invariant-mass cut can be calculated reliably, allowing for a determination 
of | Kb I at the 10% level. 



2 Construction of the hybrid expansion 

A cut q 2 > q 2 on the dilepton invariant mass squared q 2 = (pi + p u ) 2 implies 

M x <E x <^ r {M 2 B + M 2 x -ql) , M 2 <M 2 X < [m b - y^)' (1) 



for the hadronic invariant mass M x and energy E x . For g 2 = (M B — M D ) 2 both variables 
vary between M w and M D , and thus no charmed hadrons are allowed in the final state. 
If the cutoff is chosen higher (as may be required for experimental reasons), the maximal 
values of M x and E x become less than the D-meson mass. The kinematical variables 
that enter the partonic calculation of the decay rate are the parton invariant mass and 
energy, a/p 2 and v ■ p, where p = m b v — q with v being the velocity of the _B-meson. It is 
convenient to work with the dimensionless quantities £ = p 2 /{v • p) 2 and z — 2v • p/ra&, 
the phase space for which is 

< f < 1 , < z < z max = | (l - . (2) 

While £ is of order unity, independent of the value of the quantity z is parametrically 
suppressed by 

6 = — , with fi c = — . (3) 

nib 2m;, 

With the optimal value q 2 = (M B — M D ) 2 ~ (m^ — m c ) 2 the characteristic scale of the 
parton momenta is /i c ~ m c . 

The strong-interaction dynamics of the inclusive decay B — > X u I v is encoded in the 
imaginary part of the hadronic tensor T^ v defined as 

T,u( P ,v) = \ J d^e^- m ^^{B{v)\T{jl{x), J„(0)} \B(v)) 

= [ppVv + v^p u -g^v-p- ie fMUa/3 p a v 13 ) T x 

- 9^v T 2 + v^v v T 3 + (p^v u + v^p v ) T 4 + p^p v T 5 , (4) 

where the scalar functions Tj depend on p 2 and v ■ p, or equivalently, on £ and z. The 
weak current = M7 M (1 — 75)6 mediates the b — > u transition. To separate the long- 
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and short- distance physics one usually performs the operator product expansion of the 
hadronic tensor and calculates the corresponding Wilson coefficient functions in per- 
turbation theory. However, since in the present case the characteristic scale fi c of the 
momentum flow through the operator product is lower than the 6-quark mass, the per- 
turbative expansion of the Wilson coefficients involves logarithms of the ratio /z c /m&. 
To improve the quality of the perturbative results these logarithms should be summed. 
This can be achieved by absorbing the physics above the scale \x c into the couplings of a 
low-energy effective theory. 

The relevant effective theory is the HQET, and the expansion of the weak current at 



NLO reads 14, 15 



J, = Cii/i) Oi + ^Yl B ^ Qi + °i l l m D ■ ( 5 ) 

i j 

The two operators contributing at leading order are 

0\ = ulI^K , 2 = ULVfj.hu , (6) 

where Ul = u{l + 75), and h v is the effective heavy-quark spinor of the HQET. At NLO 
in 1/rrib, the operator basis contains ten independent operators. Six of these are local 
operators involving a gauge-covariant derivative on one of the quark fields: 

Qi = QLl^iPK , Q 2 = qLV^iJpK , Q 3 = qLiD^hv , 

Qi = qL{iv-D)^fJi v , Q 5 = qL^v-D^v^K , Q e = qL^D^K . (7) 

The remaining four involve time-ordered products of the leading-order currents with the 
operators Oy m = h v (iD) 2 h v and O mag = y h^a^G^h^ appearing at order l/m& in the 
HQET Lagrangian: 

Q 7 ,8 = i J d 4 x T {Oi , 2 (0), O kin (x)} , Q 9)10 = i J d 4 x T {O lj2 (0), O mag (x)} . (8) 

These nonlocal operators are present because we work with eigenstates of the lowest- 
order HQET Lagrangian, so that all dependence on the heavy-quark mass is explicit. It 
is convenient to treat the operators Q7,,..,io on the same footing as the local operators 
Qi,...,6> because the nonlocal operators mix into the local ones under renormalization. We 
have recently completed the calculation of the two-loop anomalous dimension matrix for 
the operators Qj and solved the RG equations for their coefficients Bj(fU) fL7j . 

When the expansion (|)D is inserted into the operator product in (f|), all dependence 
on the 6-quark mass becomes explicit and is contained in the Wilson coefficients Ci, Bj 
and the l/m& prefactor of the higher-dimension operators. In the next step, we perform 
the operator product expansion and obtain forward 5-meson matrix elements of local 
operators of the form (B(v)\ h v . . . h v \B(v)) in the low-energy theory. To order (A//^ c ) 2 , 
these matrix elements can be expressed in terms of the parameters [18[ 



Xl= 2M~ B {Bl ° km lB) ' XM = {Bl ° maS lB) ' (9) 
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The anomalous dimension of the kinetic operator vanishes because of reparameterization 
invariance |19] . Up to higher orders in the heavy-quark expansion, A2 is determined from 
the mass splitting \ 2 {m b ) ~ J(M|» - M|) ~ 0.12 GeV 2 . 

By dimensional analysis, the contributions from time-ordered products of higher- 
dimension operators in (H) are suppressed by powers of 6 = fi c /mb. At NLO in 5, it 
suffices to evaluate time-ordered products with operator insertions of the type {Oi,Oj} 
or {Oi,Qj}. The decay rate is then obtained by contracting the result for the hadronic 
tensor with the lepton tensor and integrating over the portion of phase space with q 2 > q 2 : 

1 Zrnax ^ ^max 

T(6)^T(B^XJu)\ q2>q2 = Jdi J <feJ^L = -L j di J dzT{£,z), (10) 

151=1 

where 2 max = 25 + <5 2 £ + 0(5 2 ) from (^[). Using the fact that the hadronic tensor is 
an analytic function of £ apart from discontinuities on the positive real axis, we have 
written the expression for the rate T(5) as a contour integral in the complex £-plane. 
The integrand of the contour integral is |20| 



r(£, z) = yfl^l \m b z [l2 - 8* - (4 - 3z)z f] T x + 6(4 - Az + z 2 i) T 2 

+ 2z 2 (1 - (T 3 + 2m fe T 4 + m 2 b T 5 ) | , (11) 

where 

^ " 1927T3 (12) 

is the tree-level decay rate. 

At NLO in RG-improved perturbation theory, the matrix elements of the hadronic 
tensor must be evaluated including corrections of order a s in the low-energy theory. 
These are obtained by evaluating the contributions of all physical cuts of the diagrams 
shown in Figure [TJ. These diagrams contain infrared divergences for £ — > from soft and 
collinear gluons, which cancel once the differential rate is integrated over phase space, 
because the rate T(S) is an infrared-safe quantity. The contour- integral representation 
in (|1(]) avoids the singularity at £ = 0. After integrating over the loop momentum, 
the function T(£, z) is given by an integral over Feynman parameters. At this point 
we exchange the order of integration and integrate over the phase space first. The z- 
integration is trivial, because the dependence of the loop integrals on z is fixed by their 
mass dimension. Note that the upper limit z max of the integration depends on £. After 
expanding to NLO in S, the ^-integrals can all be reduced to 



{G a (a);H a (a)} = ±-. /^^^{l^}, (13) 



I5l=i 
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Figure 1: One- loop graphs contributing to the hadronic tensor. The crossed 
circles represent insertions of the operators Oi or Qj. The solid squares in the 
last eight diagrams denote insertions of Okm or O mag . 



where the denominator (a— £) Q with a > is the result of the Feynman parameterization. 
Some of the diagrams have a pole at £ = due to the propagation of a single w-quark. 
The integral G a (a) evaluates to 

6(1 -a) 



G a (a) 



7T 



2r(a)r(|-a) (1 



(14) 



The second integral picks up contributions for a < 1 and a > 1. It is needed only for 
a = 1 + e and a = e (with e = 2 — d/2 the dimensional regulator). We find 



+ 



a 

6(a- 1) 



jl-v 7 ! 3 ^ l-eln4(l-a)] - 2e ln(l + v 7 ! 3 ^) + 0(e 2 ) | 



,1+e 



(15) 



We have not expanded the result for a > 1, since we will have to integrate it from 
a = 1 to oo when performing the integrals over Feynman parameters. The corresponding 
expression for the function H e (a) can be obtained by integrating the above result with 
respect to a. 

We now present the results of our calculation. We expand the rate T(S) in powers of 
5 and define 



F(5) = 16T 5 3 (f 3 + 5 f 4 + 5 2 f > 5 ) , 



(16) 
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where the dimensionless expansion coefficients T n collect terms of different order in the 
hybrid expansion. The leading-order contribution T 3 was calculated in ||; the NLO 
term r 4 is the main result of the present work. We find that 

r 3 = c 1 (n) 2 F 3 [c 1 ,c 1 ], 

f A =J2 C ^ C ^ F ^ ^ + E Cl ^) Bj(fi) F 4 [C 1 , B 3 ] , (17) 

i=l,2 j=4,6,7,9 

where {Cp = 4/3 is a color factor) 



F 3 [C 1 ,C 1 ] = l + ^M('61n^ + 19-^)-?^M + ... 

An \ 2fi c 3/2 fij 



C F a 8 (fi) ( fi 29 8vr 2 \ 1 3Ai + 3A 2 (/x) 



L ' J 4 \ 4vr V 2/x c 2 ° ; ' 1 



L ' 6J 4 \ 4vr V 2// c 2 ° ; ' 1 





1 3 A 2 (/i) 




r 4 ni 


¥)! 


I 3A 2 („) 




[ 4 /i c 2 



F 4 [Ci, -B9] = 1 Cf0iM (einJ^+ 2 4)- 9 -^ + .... (18) 



4tt V 2/i c 2/2 /i 

Note that the Wilson coefficient C 2 starts at NLO in a s , so there is no need to include 
0(a s ) corrections to the quantity F 4 [Ci,C 2 ]. The coefficients B 7 = Cj Cy m and Bg = 
C\ C mag are determined in terms of products of coefficients appearing in the leading- 
order expansion of the currents and in the HQET Lagrangian. C mag is the coefficient 
of the chro mo-magnetic operator mag , and Ckm = 1 is the coefficient of the kinetic 
operator Oki n - The remaining coefficients, _B 4 and B e , were calculated at NLO in [O. 
The Wilson coefficients, as well as the matrix elements in the low-energy theory, depend 
on the choice of the renormalization scheme. The scheme dependence cancels in the 
result for the decay rate T(5). In [|T^] we have considered a class of renormalization 



schemes parameterized by the quantity A RS = 2 ln(/ij^g//i RS ). Our results ( |18|) refer to 
the MS renormalization scheme (with anticommuting 75). The corresponding results in 
a different scheme RS are obtained by replacing 

ln /^^ ln /%s Ars_ (19) 
2^ c 2/i c 2 v ; 
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That way, it can be checked that the scheme- dependent terms proportional to Ars indeed 
cancel in the final result. 

The above expressions are given in terms of the 6-quark pole mass m& = ml° le , which 
is the expansion parameter that enters in the HQET. However, this mass parameter is 
plagued by infrared renormalon ambiguities, and it is better to eliminate it from the final 
result for the decay rate. It is well known that the convergence of the perturbative series 
of inclusive decay rates can be improved by replacing the pole mass in favor of a low-scale 
subtracted quark mass, which is obtained from the pole mass by removing a long-distance 
contribution proportional to a subtraction scale fif ~ several x A | 



24 . In 



addition to mj,, the characteristic scale \i c = (ml — q^)/{2m^) of the hybrid expansion 
and the associated parameter 5 = fi c / m b depend on the definition of the quark mass. 
It has been shown in [R| that one should choose fif = (fi c with a parameter ( = 0(1) 
so as not to upset the power-counting in the hybrid expansion. Specifically, we use the 
potential-subtracted (PS) mass m^ s introduced by Beneke [F2TJ. At NLO, the relation 
between the pole mass and the PS mass readsQ 



ni 



pole 



4cx s (//) 
3vr 



1 + 



2tt 



H f 6 



+ . . . 



(20) 



where /3 
that 



11 — \ nj is the first coefficient of the (3 function. With fif = (fi c it follows 



pole 



m 



PS 



(00 



l + 6( 



4a s (/i) 
3tt 



+ 0(a 2 s ) 



(21) 



From now on, we use the PS mass m^ 



m 



PS 



(Cfi c ) in all our equations unless indicated 



otherwise. The parameters /i c and S are defined with respect to this choice. This results 
in extra terms in the perturbative expansion proportional to (. The results in the pole 
scheme can be recovered by setting ( = 0. 

In 0, the authors have eliminated the pole mass in favor of the so-called Upsilon 



mass m^ ||22|| , which up to a small nonperturbative contribution is one half of the mass 
of the T(IS') resonance. The relation between the Upsilon mass and the pole mass is 



m 



pole 



is , 2m fe ls a,(/i) s 



9 



1 + 



a s (fi) 



71 



A) In 



mfa s ([i) 



11 



-4 



(22) 



In higher orders in this expansion the logarithms £ = ln[/i/jmj s a s (/i)] appear in the 



form [£ n /n\ + £ n /(n — 1)! + • • • + 1] ~ exp(£), which yields a factor ml s a s (jj,)]. 
Hence, schematically the relation between the masses becomes 



pole 



m 



IS 



+ /i 



6 



1 + nonlogarithmic terms] 



(23) 



which shows that despite of the unusual form of (|22]) the Upsilon mass is a low-scale 
subtracted mass in the sense described earlier. (One may, however, question whether 

1 The NNLO contribution to this relation is known, but we will not use it in the present work. 
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this argument can consistently be applied in low-order calculations.) Based on this 
observation, the authors of [^2| have suggested to treat the 0(a 2 s ) term on the right- 
hand side of ( m) as an 0(a s ) correction when the Upsilon mass is used in the context 
of perturbative expansions in the B system. Therefore, in our case we should replace 
[1 + | a s (fi) 2 + . . . ] and treat the 0(a 2 ) term as a term of order as in 



m 



pole 



m 



IS 



(pl|). This means that our formulae apply without modification to the Upsilon scheme, 
provided we set ( = 7ra s (/i)/(65) and consider this to be a parameter of order unity. 

Inserting the explicit expressions for the Wilson coefficients and matrix elements into 
(|17|), and evaluating the results for N c = 3 colors and rif = 4 light quark flavors, we 
obtain 



x- 12 / 25 r 3 
,*- 12 / 25 f 4 



l + kx 



71 

a s (m b ) 



21„iL + 4 C + fc 2 )^ ., 

2fl c J TT 2 



3 x 9/25 H m b) 



+ 



9 7T 



+ X 



-9/25 



4 fc5 « £ K) + 

9 71 



9 2/i c 



27 



7T 



+ Inx 



25 7 7T V 25 2 /i c 25 s 7 7T 



,'37 12 , 10 _ 6/25 97 
+ T7-^lni+- x 6/25 



18 25 



c 9/25 ] ^9/25 A 2(^b) + 



(24) 



where x = a s (fx) /a s (m b ), and we have used that A 2 (//) = x 9 / 25 A 2 (mb). We have ex- 
tracted an overall factor x 12 / 25 from the expressions for the quantities r„, which reflects 
the leading-order anomalous dimensions of the heavy-light currents in the time-ordered 
product in (H). The coefficients fcj summarizing the NLO corrections are 



k 7 



9403 14vr 2 



3750 225 
23161 1123vr 2 



+ 



180000 8100 
988 56vr 2 



-3.122. 



1.497, 



5625 2025 
18806 56tt 2 



15625 1875 



0.449 , 
~ -1.498, 



k 2 



23153 62tt 2 



3750 75 

22106 196vr 2 
+ 



16875 675 
661723 389tt 2 



540000 900 
45056 248vr 2 



15625 625 



-1.985, 
4.176, 
~ 3.040 
-1.033. 



(25) 



Since we have solved the RG equation for the Wilson coefficients to obtain their 
running from the 6-quark mass down to the scale \x ~ fi c , all logarithms of the ratio ///m& 
are resummed into the running couplings a s (fi) and a s (ra&). Note that the expressions 
(p4j) are scale- and scheme- independent at NLO. As a check of our result, we expand 
a s (n) in terms of a s (rrib) in order to recover the known one-loop expression for the decay 
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rate. This yields 



r. 



1 + 



Q -i(21n^ + 4 C + ^ 
7T V 2/x c 3 



8tt s 



m 6 



20 



7 



+ — -In— C + — + 

7T \3 2/x c 3 s 18 



3 A 2 
2^ 



+ . . . 



S7T 



9 



+ ... , 



(26) 



in agreement with the findings of [g, [25] . The numerical effects of the NLO resummation 
of logarithms achieved in || and in the present work are rather significant, despite the 
fact that there is limited phase space for operator evolution between the 6-quark mass and 
the scale \i c ~ 1 GeV. In the pole scheme (i.e., with £ = 0) and with q\ = (Mb — Md) 2 , 
the running from \x = down to a characteristic scale /i = /x c decreases the value of 
T3 by a factor 0.47, and it changes the value of I\ by a factor —0.49. The resummation 
effects are less pronounced in the PS scheme with ( = 1, where we find that the value 
of T 3 increases by a factor 1.10, whereas the value of T 4 decreases by a factor 0.66. 

Finally, we quote the expression for the remainder T>5 in (|16|), which contains all 
terms suppressed by at least five powers of 5, and for which we do not perform a RG 
improvement. As we will see below these contributions are numerically very small, so 
that a leading-order treatment is indeed justified. We obtain 



x 



-12/25 y 



>5 



a s (/z) 



71 



G(6,Q + {l-6) 



Ai + 15A 2 (/i) 



(27) 



where the scale \i is undetermined and will be chosen between \x and ra& in our analysis 



below. The function G(8, £) can be obtained from the results of |25| and is given by 



G(S, C) 



125 5 



L2(1-2<J) 



7T 

~6 



Vl2 



85 + 285 2 



128 



3 5 3 + 205 4 J ln(l - 25) 



1 

65 1 



+ 



4(1-^) 
3<5 2 



L 2 (25) - L 2 (l - 25) 



7T- 

"6~ 



+ 4C 



8 1 or ^ 544 76 c 

— In 25 + 4C 5 + 0(5 2 ) . 

15 225 135 v ; 



(28) 



For the optimal choice q$ = (M B — M D ) 2 , corresponding to the largest value of 5, and 
in the PS scheme with C ~ 1 an d A* c < A* < 2/i c , we find that the relative contributions 
of the three terms involving T3, and T>5 in (|I1J) are approximately 1 : —0.15 : 0.05, 
indicating a rapid convergence of the hybrid expansion, and supporting our claim that 
resummation effects are unimportant for the quantity T>5, especially since more than 
80% of the contributions to this quantity come from power corrections. 
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3 Numerical analysis 



We now turn to the numerical evaluation of our result for the decay rate T(5) and 
investigate the perturbative uncertainties in the calculation. The uncertainties related 
to higher-order power corrections ~ (A//i c ) 3 have been estimated previously |7|, ||. 

The value of the PS mass at the subtraction scale /i 2 = 2 GeV has been determined 



using a sum-rule analysis of the bb production cross section near threshold [56]. The 
result is mf s (// 2 ) = (4.59 ±0.08) GeV, which corresponds to m h (m h ) = (4.25 ±0.08) GeV 
in the MS scheme. The most recent value for the Upsilon mass as reported in the second 
paper in |27j is ml = (4.68 ±0.05) GeV and corresponds to m&(mb) = (4.16 ±0.06) GeV. 
This paper also contains a sum-rule analysis of the bb threshold region, yielding m^m^) = 
(4.17 ± 0.05) GeV. Whereas charm-quark mass effects were included in p7) they were 



neglected in |26|| . When comparing the results for the B — > X u lv decay rate obtained 
using different mass definitions, one should use numerical values that correspond to the 
same value of the MS mass. Combining the values reported above, and applying a small 
charm-quark mass correction in the case of p6f , we choose rn^im^) = (4.20 ± 0.06) GeV 
as a reference value. The correspondingly shifted values for the low-scale subtracted 
quark masses are 

ml s (2 GeV) = (4.52 ± 0.06) GeV , mf = (4.72 ± 0.06) GeV . (29) 

Below we will indicate how our results would change if different mass values were used. 

Given a value of m^ s (/z 2 ), we use (p0|) to compute the running PS mass at the scale 
jif = (fi c by solving the NLO equations 



™b (M/) = m b (/i 2 ) ± n f 



3tt 



1^2 _ -A A 203 a s (fj l2 ) \ _ 25 o; s (/x 2 ) ^ /i 2 



/Jf J \ 36 7T / 6 7T /ij 



CW=^ + V^ + i (so) 

For instance, with the optimal choice q 2 = (Mb — Mp) 2 we obtain /i c ~ 1.09 GeV, 
m PS 4.67 GeV and 5 ~ 0.23 for £ = 1. The next question that arises is which values 
of the parameter ( we should consider. A priori any value ( = 0(1) is reasonable, but 
in practice the choice is constrained by the following considerations. First, it is known 
that the characteristic scale for the total B — > X u I v decay rate, while parametrically 
of order rrib, it is numerically much smaller than the 6-quark mass, of order 1-2 GeV. 
A theoretical argument supporting this finding has been presented in Because the 



cut on the dilepton invariant mass reduces the phase space for the hadronic invariant 
mass and energy, the characteristic scale for the decay rate T(5) is lower than that for 
the total rate. Hence, the largest reasonable value of ( is about 2, corresponding to a 
subtraction scale /if = 2/i c w 2 GeV. Likewise, we should not consider renormalization 
scales larger than about 2 GeV. This conclusion is also borne out by closer inspection 
of the perturbative expansion coefficients in (|24|). Whereas for £ w 1 the coefficients of 
the various terms proportional to a s (fi) in the expressions for T 3 and T 4 take moderate 
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Figure 2: Scale dependence of the B — > X u I v branching ratio in units of \ V u b\ 2 
for q 2 — {Mb — Md) 2 , and for different mass definitions. The two bands refer 
to the PS mass subtracted at fif — \i c (dark), and to the Upsilon mass (light). 
The dashed lines refer to extreme choices of the subtraction scale for the PS 
mass. The gray rectangle shows our estimate of the perturbative uncertainty 
(see text for further explanation). 



values (and are of much smaller magnitude than in the pole scheme), some of these 
coefficients become large for £ > 2. On the other hand, ( can also not be taken too 
small. The formal limit ( — >• corresponds to the pole scheme and reintroduces the bad 
perturbative behavior of expansions based on using the pole mass. In addition, there are 
significant uncertainties in the calculation of the running PS mass when the subtraction 
scale is much below 1 GeV, where NNLO corrections become important. Based on these 
observation we will employ the PS mass subtracted at the scale fif = fi c = 1) as our 
default choice, and consider the choices ( = 2 and ( = 0.5 as extreme variations. 

To present our results for the branching ratio we define a function B{q^) such that 

Br(5 - X u lu)\ q , >ql = B{ql) x \V uh \ 2 (rj^) > ( 31 ) 

where tb is the 5-meson lifetime. In Figure we show the renormalization-scale depen- 
dence of B{q^) for the optimal choice q 2 = (Mb — Mr,) 2 and different mass definitions. 
We scan the parameters Ai and p, entering the expression for T> 5 in (|27| ) in the ranges 
—0.45 GeV < Ai < —0.15 GeV and /i < p, < m^. The sensitivity to the values of these 
parameters is very small, as indicated by the widths of the bands. The vertical gray lines 
indicate the "perturbative window" between /i = \i c and 2/i c (calculated with ( = 1). In 
this range there are no large logarithms, and the coupling a s (fi) is reasonably small. We 
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Table 1: Theoretical predictions for the branching ratio for different values of 
(?q. The optimal choice of the cutoff is e/ 2 = (M B — M D ) 2 ~ 11.6 GeV 2 . Errors 
are added in quadrature and symmetrized in the final result. 



ql [GeV 2 ] 


B(ql) 


5m b 


pert. 


(A//i c ) 3 


A(g 2 ) 


10.5 


26.9 ±4.5 


+3.3 
-3.0 


+2.1 
-3.0 


±2.0 


8.8 


(M B - M D f 


20.9 ±4.0 


+2.9 

-2.6 


+1.8 

-2.3 


±2.0 


9.9 


13.0 


14.8 ±3.4 


+2.4 
-2.2 


+1.3 
-1.5 


±2.0 


11.7 


15.0 


7.6 ±2.7 


+1.8 

-1.6 


+0.7 
-0.5 


±2.0 


16.4 



observe good stability of the results in this region, and an excellent agreement between 
the PS scheme and the Upsilon scheme for scales above 1.7 GeV. For lower scales, the 
results in the Upsilon scheme become unstable, whereas the stability interval for the 
PS scheme extends down to /x < 1 GeV. Remarkably, the curves corresponding to the 
different mass definitions all cross at a point close to the center of the window. We take 
the central value of the band corresponding to the PS mass with ( = 1, evaluated at 
the center of the window, as our central value, and use the results obtained with ( = 1, 
/i = /i c and ( = 2, /i = 2/j c to get an estimate of the perturbative uncertainty. 

The error bar at the default point reflects the sensitivity to the value of the 6-quark 
mass. The branching ratio scales as approximately the tenth power of rrib |J, and because 
of this strong dependence this is the dominant theoretical uncertainty in our prediction. 
The dependence on the 6-quark mass can be parameterized as 



ni 



PS \ A (?o) 



B ^*{l£tev) ' (32 > 

where mf s = m^ s (2GeV), and the exponent A(g 2 ) increases with q 2 (see Table [l] be- 
low). Alternatively, one may substitute the ratio mP/4.52GeV with mP/4.72GeV or 
^( m fe)/4.20 GeV. The ±60MeV error assumed in the figure is less conservative than 
the ±80MeV advocated in [23], but larger than the ±50MeV found in [27[ . 



Table [I] shows our final theoretical predictions for B(ql) as a function of the cutoff on 
the dilepton invariant mass. We separately list the errors due to the uncertainty in the 
6-quark mass, the perturbative uncertainty, and higher-order power corrections. Each 
of the quantities T n in QTBD receives corrections of order (A//i c ) 3 , whose contributions 
are essentially unknown. Because of the prefactor 5 3 = (fic/mj,) 3 in this equation, the 
corresponding contribution to the decay rate is to first order independent of g 2 and 
scales like (A/m&) 3 . Correspondingly we assign an error of ±16Po(A/mfe) 3 ~ 2 [K^) 2 , 
where we have used A = 500 MeV as a typical hadronic scale. Recently the possibility 
has been entertained that it might be possible to lower the q 2 cutoff below the charm 
threshold, because the contributions from low-lying charm states can be subtracted using 
experimental data [29]. We have therefore included a result for q^ = 10.5 GeV 2 in the 
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table. 

The results in Table p] compare well with the ones obtained in ||, where a value 
corresponding to B(q$) = 20.0 ±3.9 was obtained for q 2 = (Mb — Mo) 2 - A new element 
of the present analysis is the study of different mass definition schemes, which allows 
us to obtain a more reliable estimate of the perturbative uncertainty. The authors of 
have calculated the decay rate in the Upsilon scheme using fixed-order perturbation 
theory at the scale fi = nib. With this choice of scale and scheme we find B^q^) = 17.8, 
significantly less than the central value shown in the table. These authors have also 
included a partial calculation of higher-order corrections (terms of order fioce 2 ) and used 
ml s = Mx(i5)/2 for the Upsilon mass (with no uncertainty), thereby neglecting the 
nonperturbative contribution estimated in |p7| . Their result reported in [7| corresponds 
to B(q%) = (19.3±0.8 [pert] ±1.7 [power]). Since the Poa 2 terms are related to the choice 
of scale in the leading-order correction, it is not clear to us why the quoted perturbative 
error is smaller than the contribution of these terms. More importantly, however, the 
dominant theoretical uncertainty due to the strong sensitivity to the 6-quark mass has 
been neglected in M. 

4 Conclusions 

We have presented a calculation of the branching ratio for the inclusive decay B — > X u I v 
with a cut q 2 > q 2 > (Mb — Md) 2 on the dilepton invariant mass. Since the typical parton 
momenta ~ fi c after the cut are of order the charm-quark mass, we have performed a 
two-step expansion in the ratios [i c jm,b and A//i c and summed logarithms of the form 
a s ln(fi c /m b ) at next-to-leading order. To improve the quality of the perturbative result, 
we have eliminated the pole mass in favor of a low-scale subtracted 6-quark mass. We 
find that both the potential subtracted and the Upsilon mass lead to consistent results for 
the rate. Considering the residual renormalization-scale dependence and the variations 
between different definitions of the heavy-quark mass, we estimate the perturbative 
uncertainty to be of order 10%. This is comparable to the size of unknown higher- 
order power corrections, but smaller than the error arising from the uncertainty on the 
value of the fe-quark mass. 

If a value of the cutoff not far above the optimal choice q 2 = (Mb — Mp.) 2 can be 
achieved experimentally, we conclude that \V u i\ can be determined with a precision of 
about 10%. This would be considerably less than the current theoretical uncertainty in 
the value of this important Standard Model parameter. 
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